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Abstract 

We study the infrared divergences arising from gravitational loops in the standard cos- 
mological perturbation theory. We provide a simple solution to the problem at all orders of 
cosmological perturbation theory by redefining the perturbation theory in terms of a local 
observer. We propose to reformulate the standard perturbations in the in-in formalism, 
and obtain an infrared safe perturbation theory. Our results do not depend on any infrared 
cutoffs or similar parameters. We then present an explicit example of graviton one-loop 
corrections, and briefly discuss non-Gaussianities. 
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1 Introduction 

The primordial inflation [IJ is perhaps one of the most important paradigms of the early Universe, 
for a recent review, see [2j. Inflation is responsible for stretching the initial perturbations to 
the observable scales in the cosmic microwave background radiation (CMBR), and seeding the 
initial perturbations for the large scale structure formation [3]. Since the future observational 
constraints will provide a better understanding of the inflationary dynamics and its potential, 
it is then important to reach the desired level of accuracy by studying higher order quantum 
corrections to the cosmological perturbations. For a review on cosmological perturbation theory, 
see [1]. 

The correlators at loop level in an inflationary setup appear to be plagued by infrared (IR) 
divergences for soft quantum modes whose wavelengths are extremely large, known as the super 
Hubble fluctuations. The divergences appear when these fluctuations are summed over in the 
loops [3 El El El O [ini El [121 [13] • Regularizing the integrals via IR cutoffs does not solve the 
problem, but turns the divergences into large Logarithmic corrections depending on the cutoff 
("box size"). There has been a debate about the question if such correction are physical or not0. 

Typically, the IR divergences are a signal of an ill-posed physical question, and therefore 
their resolution depends on our understanding of the physical system and the approach we 
adopt to tackle the problem. For instance, an unphysical element in the approach which leads 
to the divergence could be due to an erroneous definition of an initial and final vacuum in a 
scattering process, as it happens in the case of soft photons or gluons emitted as a result of any 
quantum process |18j . The IR divergences may also arise if the perturbation theory has been 
wrongly organized without taking into account the relevant scales that make some contributions 
unsuppressed, as for the case of field theory at finite temperatures [T9] . 

^ There are attempts to address the IR issue by studying the pre inflationary phase in the early Universe, which 
modifies the long wavelength behavior of the solutions of the field equations, thus ameliorating or even canceling 
the IR divergences [131114) . However, this approach has a certain degree of arbitrariness, depending on the choice of 
a specific pre-de Sitter scenario. For example, a possible solution is to evolve the perturbations from a contracting 
phase to the expanding phase in a singularity free bouncing cosmology, where gravity becomes asymptotically 
free in the ultra violet regime [15]. Modification of the perturbation equations capable of ameliorating the IR 
issue (as they affect the spectral index) can arise also in models such as chain inflation, due to the interactions 
among the different components of the system [161 117| . 
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Another point of view on this problem has instead considered the dependence on the box 
size (IR cutoff) of the correlators as a physical input. The question then arises - whether the 
correlators, depending on the box size, had then to be averaged over a distribution of boxes 
partitioning the Universe on super Hubble scales, or if the correct physical interpretation would 
be to fix the box size to the desired scale of observation and keep the Log-enhancements H. 
Both these approaches have negative aspects - the first one does not take account of the fact 
that the observer is not capable of averaging over boxes larger than its Hubble patch, while the 
second approach still suffers from the ambiguity on how to define the box size (for example, 
via comoving coordinates or physical ones), which therefore makes the result for the correlators 
ambiguous. 

In this paper, we will point out that the issue of IR divergences in gravitational loops in 
the cosmological correlators can be solved by defining a local observer who is responsible for 
measuring the observable quantities. We will argue that once the observer and observables are 
well defined, the IR divergences will turn out to be an artifact of some unphysical assumptions, 
usually taken in the definition of the perturbation theory. 

Our implementation of the principle of locality will be different from what discussed in other 
studies in the literature El- In fact, our resolution is to redefine the standard perturbation theory 
in the in-in formalism based on one physical principle: 

• any observable quantity should be defined in terms of a local observer who is measuring 
those observables. 

As a result, the observable quantities are shown not to have IR divergences, nor to depend on 
the IR cutoffs. 

The paper is organized as follows; in section [21 we start with a brief introduction of the 
standard in-in formalism and the quantization of the perturbations. We then provide a general 
discussion of the IR behavior of loop correlators in section [31 using the example of self-energy 
diagrams, and discuss the ambiguities in the cutoff procedure for regularizing the integrals. We 
present our resolution of the IR issue in section [H First, in section I^TTl we discuss the nature of 
the true observables and how the in-in formalism, as it stands, fails to account for them. Then, in 
14.21 we show that all IR divergences at all orders arise only from certain specific contributions to 
the correlators, and can be traced back to the definition of the correlators using the background 
unperturbed metric; finally we discuss how to eliminate the IR divergences by amending the 
standard in-in computations in section 14.31 and provide a detailed example at one-loop for a 
specific case in section 14.41 Finally, we will briefiy discuss the non-Gaussianities in section 14.51 
and conclude in section \5\ 

2 Basic formalism 

The appearance of large IR corrections in the perturbation theory in an inflationary background 
is directly linked to the presence of super-Hubble fluctuations of light flelds [2il [25l[i6l 12711281 129] . 
it is therefore not related to any speciflc property of the fleld, i.e. graviton or scalar field, running 
inside the loops. Here, we will briefly introduce the basics of perturbation theory as it is usually 

^ This also seems to be the point in Ref. [201121) . as the authors distinguish between IR and non-IR perturbations 
with respect to some typical observer's scales qo and L^^ , and therefore the spectrum which they claim to be 
IR-safe would then depend on this cut-off procedure. We comment more on this in section 231 

^The issue of locality has been discussed, with a different approach, also in [221 123) . Their proposal, however, 
entails considering non-standard gauge transformations that are singular at large scales and results in new gauge 
conditions [52] or in the necessity of using perturbation fields different than the usual curvature and tensor ones, 
with a prescribed set of boundary conditions, in order to avoid the IR divergences [23) . 
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formulated using the in-in formalism and the path integral quantization, for a review, see 0; at 
the same time, we list our conventions. 

The cosmological perturbation theory begins with the definition of a threading and a slicing 
of a spacetime (foliation) in the unperturbed background metric, which is fully homogeneous and 
isotropic. Then all the fields (including the metric) are written distinguishing the background 
and the perturbations according to the chosen threading and slicing, for a review, see The 
background is considered to be classical, while the perturbation fields are quantized. 

Let us consider a scalar field ^> (similar procedure will follow for tensor fields, gravitons, 
once a polarization tensor is specified), and write it in terms of the background (poit) and the 
perturbations (j){x, t) as: 

^{x,t) = Mt) + <P{x,t). (1) 

The field (pix, t) is expanded on a basis of eigenfunctions Yj^{x) of the Laplace-Beltrami operator 
with eigenvalues —k'^. Our conventions are: k is the comoving wavenumber related to the 
physical momentum p by 

k = ap, (2) 
and the conformal time is defined via the background metric 

ds^ = -df + a^Sijdx'dx^ . (3) 

as 

dr] = a'^ dt , (4) 

where a is the scale factor. We set the reduced Planck mass, Mp = 1. The field can then be 
promoted to a quantum field, written as0: 

4^{x,r])= / : — nr Yr{x)uk{r])Zr + Yi{x)ul{'q)a\ , (5) 



(2vr)^ 



and quantized in the usual way, 

[4, %,] = <5(fc-fc')- (6) 
The Whightman function is defined as: 

H^(x,ry,x',r?') = {^\^{x,r^)^{x' ,r,'m = j ^Y^{x)Y^ {x)uk{v)uUv') , (7) 

where \Q) is the vacuum state. For convenience, we adopt planar coordinates, i.e. Yj^{x) = e^^'^, 
and, as a shorthand notation, we call Whightman function also 

Wkiv,v')=Mv)nUv')- (8) 
The power spectrum, P^, is related to the Whightman function. In particular, at late times 

P, = ^liuiWk{v,v)- (9) 

The mode functions which are entering the expansion of the field are such that Uk{r]) = Hk/z, 
where /i^ is the solution of the field equation: 



^fc +{k^--]fik = 0. (10) 



z" 



*We will denote quantum fields, both fundamental and composite, with' 
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Figure 1: Correction to the scalar two point function from an intermediate scalar and graviton. 




Figure 2: Correction to the scalar two point function from a graviton bubble. 

The quantity, z" jz^ depends on the background, it is often called the gravitational source term 
for the seed perturbations [4]. 

In a (quasi) de Sitter case, at leading order in the slow-roll parameters, the general solution 
of the mode equation is given by: 



uM = ^ = C^HW^-krj) + c^Hl'H-kr,) . (11) 
z 2z 2z 

where H^^' are the Hankel functions. The flat space solution (Bunch-Davies vacuum) is re- 
covered for T] — )• — oo, if ci = e^^'^^a)^, C2 = 0. 

In the in-in formalism, the expectation value of any operator O at a given time rjQ is given 

by: 

{n\d{'no)\^) = {n\f (^e^i'--'^''^^) d{rio)T (^e-*^'--'^''^^) \n) . (12) 
In the case of the two-point function, expanding in the interaction Hamiltonian leads to: 

A{ki,k2,'n) = Atree{kl,k2,r]) + Aloops{kl , k2 , T]) , (13) 

where the correlators, A, depend on kf and indices are contracted using the background metric. 
The time r] is also defined by the unperturbed metric. 

3 Infrared divergences 

We discuss here the general features of IR divergences using the example of the lowest-order 
gravitational loop diagrams that we are typically interested in: the self energy diagrams, such 
as those given by the cubic interaction vertices, and the "bubble" diagram. We present these 
two diagrams in Figs. ([T]) and ([2]). The fields running in the loop can be different, however a 
generic theory with gravitation will involve a graviton loop. 

The computation of these diagrams has always been performed in a pure de Sitter space (SJ 
El El [HI 121 [lOl [HI [121 US]- However, some authors have also discussed the case of a quasi de 
Sitter background [TH [13]. In order to discuss the IR divergences, we will not need to know the 
whole expression for these diagrams. Those can be found, for example in Refs. [SI [9l [lOl \TT\ . 
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Furthermore, our observations are largely independent of de Sitter or quasi de Sitter. Let us 
briefly list here the main differences. In a pure de Sitter case the scalar metric perturbations are 
pure gauge, and the only physical perturbations are the tensor ones, i.e. that of the gravitons. 
In Eqs. (fTOl [TT]) . u becomes equal to 3/2 in a pure de Sitter case. 

In a quasi de Sitter case the de Sitter scale invariance is broken, which results in a modifica- 
tion of the term z" /z in Eq. (jlOp . It is evident from the experiments that the breaking is very 
small, within the observable ranges of scales relevant for the CMBR measurements, which corre- 
sponds to roughly 7 e-foldings of inflation. It can be parametrized with the slow-roll parameter: 
e = —H/H^ . The solution to the field Eqs. (jlUI [TT]) now has u = 3/2 + (1 — n)/2, where n is 
the spectral index [^j. For a red-tilded spectrum, as the one observed for scalar perturbations, 
i.e. n — 1 < 0, the IR divergences of the loop integrals become worse than that of the pure de 
Sitter case. 

The computation of the diagrams is quite involved. Their structure is not immediately 
transparent for what concerns the physical interpretation of the IR divergences. In particular, the 
diagrams have two types of integrations - one over the (conformal) times for each of the vertices, 
and the other over the loop momentum. There are two possible sources of IR "divergences". 
One comes from the time integrals and it is present only for some kind of interactions [7J . In our 
scenario it will not appear, while, in theories where it does, it is believed to be cured in realistic 
models by using the dynamical renormalization group techniques |10| (see also [7]). The other 
IR divergence comes from the momentum integral and it is always present. We will concentrate 
on this latter one. 

We briefly discuss the issues regarding the regularization ("box approach") of the IR diver- 
gent integrals. This analysis is actually important to understand that when the dependence on 
the cutoffs is not eliminated, the IR issue is not fully understood, although the corrections can 
be made sufficiently small. The important point here is that there is no unique choice for the 
IR cutoffs, rendering the results for the observable quantities ambiguous. 

For instance, the main difference for the IR divergence in the momentum integration arises 
from the choice of imposing a cutoff on the physical or on the comoving momentum. To under- 
stand its consequences, we consider the example of a pure de Sitter case and the simplest IR 
divergent integral, that for example comes from Fig. [2) 

A('^)(r?) = ^ / d'kWk{v,v)- (14) 

By choosing the cutoffs on the physical momentum, Puv/IR — ^UV/lR-> where Mjjy/jj^ does not 
depend on time, as in [10], one finds that for small ij 

A/ ^ 1 r^^^dk. H\ (Muv\ ,,,, 

The rationale behind this choice is that this kind of cutoffs does not break de Sitter scale 
invariance. 

Instead, choosing the cutoffs on the comoving wavenumber as, kuv/iR = ^uv/iRi where 
Km is independent of time, as in Ref. [H], and Kjjv = a{r])i2, where is the renormalization 
scale, one finds 

^Which is a function of e and other small scale breaking parameters. The precise values of the spectral index 
for gravitons and scalars are different, see [2]. 
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where is the number of efoldings from the beginning of inflation up to time r] and we have 
neglected the UV contribution proportional to log jj. This kind of cutoffs does break de Sitter 
scale invariance. The rationale behind this choice is that we integrate over modes that were 
sub-Hubble at the beginning of inflation, Km = aiHi. 

We see that Eqs. (fT5t [T6|) give very different results. First principles do not instruct us to 
prefer a cutoff either on the physical or the comoving momentslj, and thus the scale dependence 
of the "observable" changes depending on the cutoff, which clearly makes it unphysical. The 
fact that the correlators grow with the IR cutoff, that is the "size of the box" L = KJ^ also 
creates difhculty in approximating the ensemble averages via the spatial averages, which are 
performed in practical observations. In fact, the RMS deviation between the spatial and the 
ensemble averages goes like APfc ~ Pfc/(A;L)3/2 [30J, in the case of the spectrum, and, since 
grows with L due to the IR divergences of the quantum loops, it could become non-negligible 
and should be taken into account for precise measurements and predictions, for not too large 
boxes and for certain scales (in particular for scales not too different from the "box size"). The 
solution to this problem is to apply the ergodic theorem only to the true IR-safe correlators. 

To resolve the issue of IR divergences, we must then propose a recipe for the definition of a 
sensible perturbation theory that does not have any of the ambiguities we have presented. 

4 Eliminating IR divergences 

Our solution of the infrared issue will be centered on the concept of local observer, as we said in 
the introduction. In particular, our implementation of the principle of locality, will be different 
from the one dealt with, for example, in [22^ i23j. which consists in using gauge transformations 
that are singular at large scales. It will also differ, as we will see, from the approach in j20| I21j. 
which, in particular, uses an explicit cutoff q^, of the size of the typical scale of observations, 
when proposing how to resum the contributions from scales larger than q^^ 0. 

Our presentation is divided in three parts: in section 14.11 we investigate the physical point 
at the origin of the infrared divergences, by discussing what the true observables are and why 
the in-in formalism, as it stands now, does not properly account for them. In section 14.21 we 
demonstrate that all infrared divergences from momentum integrals in gravitational loops have 
a common origin, which is precisely the unphysical element of the formalism presented in section 
14.11 Finally, in section H31 we outline our solution of the infrared issue proposing a more physical 
redefinition of the in-in formalism. A detailed example of how our proposal is implemented is 
given in section 

Our notation will be as follows: we will write the background FRW metric as g^^: 

ds^ = g^^ydx^dx" = -dt^ + a^Sijdx'dxK (17) 

and, in our conventions, the background quantities, which are contracted using the background 
metric, have no labels: for example k"^ = 6^^kikj. 

The perturbed quantities are instead indicated by "over-bar". Thus, the local classical metric 
is given by (neglecting the vector part): 

= g^^dx^'dx" = -N^dt'^ + a{tfhij[dx' + lTdt){dx^ + Wdt) (18) 

^There are different valid physical motivations suggesting different kind of cutoffs, for example finiteness of 
inflation for the comoving one, or causality and/or superhorizon scale invariance for the physical one, but there 
is not an utterly univocal principle that truly stands above the others. 

'^As we will discuss at the end of section 14.31 this proposal is therefore an improved version of the "physical 
box size" cutoff technique that we discussed in the introduction: effectively the result is equivalent to just using 
from the start an infrared cutoff equal to go, with all the conceptual problems we outlined. 
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where N^,N are the shift and lapse functions, and 



(19) 



The quantities contracted with this metric will have an "over-bar", for example k = h kikj. 
4.1 True observables and in-in formalism 

We will here argue that the usual formalism for the perturbation theory does not define correctly 
the observables from the point of view of a local observer, and this is the origin of the IR 
divergences. 

Let us look closely at the definition of observables in the in-in formalism. For definiteness, 
let us consider the two-point function as a working example (our considerations apply also to 
A^-point functions, as we will discuss in section r4.5|) . The definition of the observables in the 
perturbative in-in formalism is given by the right-hand side of Eq. ()13p . The important element 
in this formula is that the observables are defined in terms of the background quantities, in 
particular the background metric (entering for example in k'^). However, the background metric 
has no physical meaning and we claim that this is the reason of an appearance of the IR 
divergences. 

In fact, any local observer uses true local clocks and ruler. The observer has no notion of 
a "background metric" and "perturbations" on it, but instead he/she uses the local metric for 
his/her measurements. We will show that when writing the perturbative correlators from the 
point of view of the local observer, that is using the local metric, the IR divergences are cured. 

In more details, the IR divergences arise for the following reason. In the in-in formalism, 
the perturbations (Ct) are promoted to quantum fields (C)7)) and path-integrated using the 
background metric of Eq. ()17p to contract the indices. We will show that in the correlators at 
loop level, all IR divergences arise when at least one of the momenta running in the loops lines 
goes to zero. In that limit the correlator becomes effectively disconnected and leads to the IR 
divergences. We represent this, schematically, in Fig. [3l 

We claim that it is precisely these operations of expansion over the background metric and 
quantum averaging of the perturbations in effectively disconnected contributions that make no 
sense from the point of view of an observer. These are the unphysical elements of the standard 
in-in perturbation theory, because a local observer will not be sensitive to the background metric 
Eq. (fT7|) , but to the local metric Eq. (fTSj) , whose classical value in the quantum theory is 



We will show in section 14.31 that by correctly defining the correlators in terms of the local 
classical metric Eqs. (jlSp . ()2ip these disconnected pieces will be removed and the infrared issue 
cured. Our IR-safe correlators will be different from those defined in other proposals, as for 
example in pO t l2T] or [22 | [23] . 

*The shift and lapse functions are Lagrangian multipliers, determined by the constraints as functions of the 
perturbations (",7 defined in Eq. (|19p . In Eqs. (|18p . (|19p . we made a gauge choice, for what concerns the 
parametrization of the metric and the choice of coordinates (x,t). In particular, we distinguish scalar and true 
tensor part, so that 7^ = d''jij = 0. The choice of a gauge does not introduce any element of non-physicality or 
arbitrariness in the description of the observable. 




(20) 



where, in particular 




(21) 
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f fd\d^KW„FK\W\5{q + K) -)• fd\n],^o fd:'KFK\W\S(K) 



Figure 3: A schematic representation of the loop corrections in an IR Hmit as shown in the text. 
A generic correlator can be written, taking care of the time-ordering, as a product of Whigtman 
functions integrated over the momenta running in the internal lines, with delta functions at 
vertices. In the IR limit of a diagram line (for example, here in the picture, a graviton's one), 
the total correlator factorizes into the product of a disconnected two-point function at small 
momentum times the rest of the correlator with a certain "vertex" . The disconnected two-point 
function (when (7—7-0) leads to the IR divergences. We will show that when redefining the 
correlators using the local metric, the disconnected contributions disappear. 



4.2 IR divergences at all order 



We will now prove that indeed all IR divergences from momenta integrals originate in the 
unphysical expansion of the local metric in perturbations and their path integration in the IR 
limit where higher order correlators become disconnected, as we have claimed in the previous 
section. Here, we will focus on the two-point function for scalar perturbations, but our analysis 
can be extended to the general case of A^-point functions and to gravitons. 

Let us first review in some detail the path-integral formulation of the in-in perturbation 
theory, for a general reference see [6]. The Lagrangian we are keen to discuss is the one of a 
nearly massless scalar field, a, which could be the inflaton, minimally coupled to gravity. 



S=ll V-det(5) [R - 9^,d^<yd,a - 2V{a)\ , 



(22) 



where V{a) is very flat, otherwise arbitrary. Using the ADM formalism | |31j . the action becomes 



S 



_det(7^) y 



5(3) 



NB}-'' - 2NV{a) + N ^ {e\E' j - {E\Y ) + 



where 



Eij = - (h 



V,;A^, 



(23) 



Vj is the three-dimensional covariant derivative calculated with the three- metric hij] and R 
is the curvature scalar calculated with this three- metric: 



-d(3) 



-2 -2C 

a e ^ 
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All spatial indices i, j, etc. are lowered and raised with the metric hij and its reciprocal. In 

particular, by computing the derivatives, such as hy in Eq. ()23p . one can find a kinetic term 
for C and 7, where indices are contracted by hij. We choose the gauge where the scalar field is 
homogeneous, and C 7^ 0- 

The two-point function in the standard in-in perturbation theory (using the path-integral 
formulation) is obtained by expanding the metric hij in terms of the perturbations (^,7. In this 
section we will not indicate quantum fields with ^ in order not to clutter the equations. The 
formalism also requires doubling the field degrees of freedom in order to account for the time 
ordering and anti-ordering in Eq. <\12\i : (" — )• {(^+,^_}, and similarly for 7. It will appear to be 
convenient to change basis to Cc = ^(C+ + C-)? Ca = C+ ~ C-i and the analogous for the 7. 

The perturbative expansion in terms of the perturbations C, 7 is obtained by functionally 
Taylor-expand the exponential e* '^''in a:\/ dethC{Cc,ic,CA,iA) ^.j^^ path integral in powers 
of the perturbations (c, Ic-, Ca, 7a: 

1^1 ^^^'"'^'^*'""^'/^^C^CA2?7c^^7ACc(^i,r?)Cc7(x2,r?) 



n,m, a=0 
u,v>0 



fe=0 



n(/ ^^^(c)^^(c)Ca(^(c),?7(c))- 
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c=0 

V 



d=0 



57a {x(d),V(d)) 



jn. dT j (fxy/ -deth C(C,c nc ,CA,lA,hc ,hA) ^ / / x 



where C is the closed time path, defined by C+l'?!!!) = C-(^in) (similarly for other fields), and the 
vacuum function is ^'vac,c/A- Finally, >C(Cc, 7c, Ca, 7A, /^c^a) = £(C+,7+>+) -'C(C-,7-,^-), 
and, after the functional derivations, Cc>7CiCa>7A — ^ 0, hij — )• in C for the usual Taylor 
expansion around the unperturbed background. 

We are now ready to investigate what kind of IR divergences from momentum integrals are 
present in the cosmological perturbation theory in a (quasi) de Sitter background. 

We find that all these IR divergences, at all orders, appear when a momenta running in a 
loop line goes to zero, as we claimed. The proof of this is as follows: in the perturbation theory 
over (quasi) de Sitter, the Whightman functions go as E|: 



1 



1^^ \k-^^ (|fc|2 + |gl2_2|fc||glcos0)i ' 



(25) 



and therefore this would diverge only when p = 0. This happens when the momentum running 
in the line goes to zero. 



More precisely, in the (quasi) de Sitter case, they would go as ~ IpI -(*-"', where n is the spectral index. 
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In particular, there are no coUinear divergences arising when two momenta are parallel, 
because the denominator of the Whightman functions (from which all other propagators can be 
obtained) never goes to zero in that case. This is very different from what happens in a scattering 
amplitude on Minkowski background, where such collinear divergences are present, for example 
when a soft gluon (or graviton) line originates from a massless line yielding a contribution to 
the diagram of the form: ~ = i^^_^yi^f^_^-^ = _2fcMg ' which diverges when A;^ is parallel to 

The IR divergences in momentum loops in perturbation theory can therefore all be recovered 
by taking the small momentum limit of the propagators in the relevant loop lines. Now we 
want to prove the other claim of ours: that all the IR divergences appear from disconnected 
contributions to the correlators. The relevant propagators in the usual ± basis are: 

G^^{r], r]') = W{x, ij, x' , r]') , G^^{r], r/') = W{x, r]' , x' , r]) , 
G++(x, y) = 9(7] - r?')G-+(r?, tj') + 0(r?' - r/)G+- (r?, r?') , (26) 
G—{x, y) = 0(7? - rj')G+-{rj, rj') + 9{rj' - v)G-+{7i, i]') , 

where W{x,r], x' ,r]') is the Whightman function (recall Eq. ([7])). Note that G~^~^ + G = 
G~^ + G ^. In the new, more convenient, (C, A) basis the correlation functions are 



The advanced and retarded propagators are related by Gyi(x,y) = G^^y^x) and vanish in the 
coincidence limit. The convenience in using this basis descends from the fact that in the IR 
limit (small momentum), the propagators behave as follows 

Gc{q^\g],r^,r^') W{q,r^,r^%^o ^ , Gn{q,r^,rj') ^e{r^-rj') ^' / (28) 

and therefore the IR divergences coming from the vanishing of the momentum in loop lines are 
accounted for by the Gc propagator only. We concentrate therefore on the expansion in CciTc- 
Let us then see how all these IR divergences arise from disconnected contributions to the 
correlators. It is easy to do it now that we have shown that they all arise when one (or more) 
momenta in a diagram line goes to zero. Indeed, the small momentum limit of the propagators 
can be found using the expansion of the perturbation fields given by Eq. ([5]), where - we recall 
- we have chosen planar coordinates: ^^(2^) = e^^'^. In this limit, i^g(2^) ~ 1, and thus the 
momenta of the lines sent to this IR limit drop out of the delta functions at the interaction 
vertices in Eq. (j24p . Therefore, in the IR limit the contributions of these lines to the two-point 
function become disconnected, and Eq. ()24p reads, once written in momentum space, as: 



IR, all orders" 



n,m7^{0,0) ' '0<a<n ^ a=0 6=0 

0<6<m, 



(IR) 



VCc'DCAD^c'DjACcikuv)Ccik2,r]) fl( [ d 



nl I ^ \ jn. dr J d^x\/-dethC{Cc ,7C ,Ca ,7A ,hc ,hA),j, / \,t, / n 

/ " X(^h)- r e "m Vvac,c(^inj^vac,A(^/inj , 

(29) 



b=0 
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where again, having Taylor expanded, CC)7C)CA)7A 0, hij — )• 5ij, in C after the functional 
derivations. The label ^jj^^ indicates the small momentum limit of the propagators (infrared 
divergent). 

We have thus proven that the equation ()29p comprises all the IR-divergent parts of the 
two-point correlator from loop momenta integrations, and these appear as disconnected contri- 
butions. With a more schematic and shortcut notation, we can write Eq. (j29p as 

{Cikl,1l)C{k2,r]))m, all orders = 

E ;^(nC):a(n7H.,.(.,):K^^3^^ ^)C(/=2, r?)) . (30) 

n,m^(0,0) ■ ■ «=0 v=0 H '-h^{v)j{v) 

With the suffix ^ we indicate that the derivatives do not act on the external fields, C{ki,r]), ({k2,r]), 
of the correlator, but on the fields entering the correlator via its dependence on hij = e^^{e^)ij. 

The form of Eq. ()30p makes it even more evident that the IR divergences are a result of 
the expansion of the local perturbed metric over the long-wavelength perturbations and the 
path integration in the limit where higher order correlators become disconnected. This result 
completes and extends that in Ref. [11], where it was shown that some IR divergences in cer- 
tain one-loop in-in computations could be recovered from semiclassical expansions around the 
perturbed metric 0. We have here proven, using a rigorous field theory formalism, that all the 
IR divergences in momentum integrals in gravitational loops at all orders are accounted for by 
the equation ()30p . The demonstration applies to A^-point function for all A^'s as well. 

Since now we see that all IR divergences in gravitational loops have the same origin, whose 
nature we can understand, we will propose a procedure to fully resolve them on the basis of the 
physical concepts discussed in section HTTl In the next section we turn to this point. 



4.3 Defining a different in-in perturbation theory 

We set out now to define a new perturbative expansion in the in-in formalism, capable of curing 
the infrared issue. We will outline our recipe, and comment on how it differs from other proposals 
concerning the IR issue. For definiteness, we will be illustrating our proposal on the two-point 
correlator {({ki,r])({k2,ri)), but, just as before, our considerations also apply to (77), or any 
other correlator of N fields, for all A^'s. 

In the usual perturbation theory around the background metric, the two-point correlator is 
given by a series 

{Cik,,r])C{k2,v))=A{k,v)Sih + k2) = {A{k,ij\^^^+A{k,ij\^^j5ih+k2) (31) 

We stress that we do not take any infrared limit of sorts here, but instead we are considering the 
full perturbative series {A{k,r])^^^^^ = A{k,ri)-^ ^^^^ + -^(^j ^)2-ioop + • • •) with the loops contribu- 
tions including the ultraviolet, finite and infrared parts altogether (as we are not concerned here 
with the ultraviolet divergences, we will assume that all correlators here and in the following 
are suitably renormalized). 

The background metric enters in this formula because it is used to contract the comoving 
momenta as (recall that we use the ADM formalism) 

~ ^blkgr ^j^j ~ ^j^j 

where we have used h^^ = 

bckgr 

^''in Ref. |21) . it was also argued that some loop IR-divergence could be obtained in this way, but the 5N 
formalism was used instead, which was questioned in 
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Motivated by our results and understanding in the previous sections, we set out to define 
a new perturbation theory, where scales are not defined by contracting with the background 
metric, but with the field 

7^,^. = (e2C(-'*)(e^(^^'*)),,). (33) 

The idea, let us repeat it once again, is that in this way we should be able to take into account 
the actual local metric that defines scales for the observer. We want to see if in this way the 
infrared issue will be cured or not. Please, observe that in this definition, the notation ( ) 
indicates a quantum expectation value in the in-vacuum and that C, 7 are the full quantum 
fields. Therefore the definition of h is the standard definition of the classical field as the in- 
vacuum quantum average of a quantum operator. In particular, it is not based on any infrared 
limit or large scale averaging or any infrared/large scale definition/quantity/procedure. 

As the background metric was entering the correlators by contracting the wavenumbers k, 
the new field h (the local metric) will have to be used to contract the wavenumbers in place of 
the background metric, as 

t = (hf kjkj = (e-2C(^'*)(e-^(^'*))^^') hkj. (34) 
Note that this quantity is very different from the one defined in Refs [21] , which is 

ti{qo,Lf = e-2^-(^'*)(e-^-(^'*))yFF ; (35) 
the latter quantity is in fact defined by using the large scale behavior 

Q,r/l^r = (2vr)-t j d^ge^^^^XfcVTfc , (36) 

where go is a suitable cutoff defining an infrared/large scale limit where the fields become clas- 
sicajii]- As we said, we use instead the full quantum fields and the usual quantum vacuum 
expectation value, see equation ([33]) . The logic at the basis of the proposal in Refs [201 [21] and 
of equation ()35p is the definition of a different background, characterized by the scale go, taking 
into account superhorizon modes, on which to do a new perturbation theory. Their power spec- 
trum is defined as the Fourier transform of a spatial average at two points xi and X2 = xi + y 
with q^^ ~ y and then distances are rescaled with the new background metric, which is not fully 
local as well (it depends on q^^ ~ \x2 — xi\). 

We, instead, put the physical concept of locality at the basis of our redefinition of the 
perturbation theory, asking what a local observer would actually use: (j33|) is not a background 
metric, but the local metric that the observer would use to define scales and a local quantities. 
Further differences among the other proposals and ours will soon appear. 

At first sight, at this point we can construct two kind of correlators using h that we might 
think to use to define a new perturbation theory: 

A!{k, 77) = M{k, 7y),_,_, + M{k, rj\^^^^ , (37) 

or 

A(k, rj) = A(k, r/),^^^.,_, + A(k, r?),„„^^ . (38) 

Here, the first object A'{k,r]) is the transformed two-point correlator under the transforma- 
tion k ^ k acting as a change of coordinates. The second object A{k, rj) is instead the non- 
transformed correlator evaluated on the newly contracted k. We claim that the latter is the one 
we should use for the new perturbation theory. 

^^In Refs [201 121) . go is the scale of observation. There is a similar quantity -defined however in position space- 
in [23] where qo is given by the Hubble rate. 
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In fact, it is straightforward to see that ()37p does not give any new perturbation theory, 
because A'{k,rj) = A{k,rj), since the correlators transform as scalar under the transformation 
k ^ k, and thus we would be just re-writing in new coordinates the same old perturbation 
theory, with all the same problems. 

Instead A{k, rj) is obviously different from A{k, rj) and thus defines a new perturbation theory. 
It is straightforward to obtain the explicit expression of A{k,r]) in terms of A{k,rj), by using 
standard techniques in calculus and the action of diffeomorphisms. We briefly illustrate the 
passages: 

• we take the correlator A{k, rj) (which comprises the tree-level and all loop contributions, 
see equation ([3T]) ) and transform its fc-dependence into a A;(A;)-dependence using 

A;2 = (h-^y^ k-k- = (e2C(^'*) (e^(^'*))^i) kikj, (39) 

• we then expand in Taylor series: 

A{k(k),rj)=A(k,r,)\ + -]^((0")((T^r)9a9-^(^,ry)|c^ 0, (40) 

n,m>0 ik^=k ) 

{n,m}yt{0,0} ^ ' 



we finally note that the first term on the right hand side of Eq. (|4U|) . i.e. 

A(k,v)\^_^_, = A{k(kU) - i((c)")((?^r>5c;5- ^(A:,,,„^=,^ 

'h r„ ll'.IILl 'h ran 



n,m>0 ■• ■■- " ip,2_p,^\ 

is precisely the non-transformed correlator evaluated on fc, which we were looking for. 
Once again, we stress that we have not taken any infrared limit or similar, and that all 
higher order contribution include all ultraviolet, finite and infrared parts. 

In fact, at this point it is not yet clear that the infrared issue has been cured in this way. In 
order to understand that, we need to look more carefully at the loops contribution in A{k^ rf) 
to assess better how the new formalism, centered around /i, differs from the old one centered 
around the background metric. Thus, by substituting ([3T]) . ([37|) in equation dH]), we find that El 

AK^\....=A'{KV\..,^ - Y i((Cr)((7^^r)555-^(A:,r?)|c,_^ 0. (42) 

fill II. ' h __9 

n,m>0 (i^2_p,^\ 

{n,m}^{0,0} ^ ' 

Clearly, A{k, rJ)^^^^^ is different from zero, because the first term on the right hand side is not 
equal to the the second term (the sum over disconnected contributions), and thus their difference 
is not null. In fact the first term is the standard loops contribution, not given by disconnected 
diagrams. Recall that all terms at each order include their ultraviolet, finite and infrared parts 
altogether. 

However, because of our results in section 14.21 we see that the infrared divergent parts do 
cancel between the first and the second term on the right hand side of equation (|42|! [^ and thus 



^Here we have also used A'{k,ri) = A{k,ri) to write everything more neatly in terms of k. 

^Note the crucial minus sign in front of all the higher order terms in our equation (|41|l . Such minus has nothing 



to do with an expansion in fc of fe^ = (e 2?(^'*)(e ^(^•'))'-'^ kikj, which would instead give 

^'(fc(fc),^) = Afc,^)|c,=7,=o+ E ^((0")((7^r)Ca".-4(fc.'7)ic,=7,=o. (43) 

because Wick's theorem always forces n and m to be even in order to have non-zero correlators. This also shows 
that (|41|) is not equivalent to a gauge/coordinate transformation. 
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A{k, rf) is free of infrared divergences: looking at equation ()30p it is straightforward to see that 
infrared divergences/Log-enhancements, and only those, do cancel out exactly order by order. 

We stress that, instead, the ultraviolet (renormalized) and finite terms do not cancel between 
the first and second term and so the loop contributions A{k, in the new perturbation theory 

around h are 1) different from those in the standard perturbation theory using the background 
metric (thus the new perturbation theory is a non-trivial modification and in principle testable), 
and 2) have the bonus that the IR divergences are absent (without using any IR cutoff procedure 
to define the new theory). 

As a comment, let us observe that, being fully expressed as a combination of usual standard 
coordinate-transformed correlators, as visible in equation ()4ip . also the new correlators will be 
sharing the same properties of the old ones (such as conservation and alike) at large scales (large 
k). Furthemore, since the field C is fully gauge invariant and, at the level of field equations, 
(j34p acts as a coordinate transformation, will be conserved on superhorizon /c-scales. Note 
also that we could have obtained our result using different gauge-fixed quantities (such as the 
inflaton perturbation) or even gauge-invariant variables such as the Bardeen potentials - in all 
cases the physical interpretation of the IR-divergences linked to the metric/clocks of the local 
observer is the same (a local dressing and rescaling of the metric or Hubble rate, see section 
SI and the prescription for the IR-safe correlator, Eq. (j4ip . changes only in the fact that the 
chosen perturbation variables must be used in place of C, 7- Note that the action is indeed gauge 
independent, and therefore the vertices are well-defined. 

As a final comment, it can be useful to compare our proposal with the one in Refs |201l21j. As 
we have already shown that their proposal is based on a quantity n{qQ,Vf' = e"^^*'' {e~'^*'')ijk^k^ , 
see Eq. ([35]) . defined using qq, L to define the classical fields Qr, lir via the IR or large scale limit 
given by Eq. ([55p . We do not use any large scale limit in our case, but employ the quantum 
expectation value of the full quantum fields in our Eqs. (p3]l . (pl|) . 

Moreover, the IR-safe two-point correlator is also defined in a different way than in our 
case. Refs. |20^ [21] define it via a transformation k — )■ K,{q, L) of the Fourier Transform of a 
spatial average at pair of points: 

A'{K{qo,L),7])tree- (44) 

In other words, their proposal is to fix a certain physical cutoff qo corresponding to the typical 
scale of the observation ( qQ^ ~ |x2 — with X2,xi the pair of points of the average), and define 
as "infrared" all wave numbers k < qo, essentially decomposing the fields as ( = Cii^'^° + C'^^'^", 
promoting only (^^>io to quantum fields (analogously is done for 7) and declaring the infrared 
safe correlators to be just the tree-level result written in the new variables, see Eq. ()35p . 

The first question that arises here is: since the correlator in Eq. (j44p is only the tree-level 
contribution, supplemented with the IR effects in Eq. ([35]) . what would be the IR-safe loop 
corrections to it in the proposal of [201 EI]? The paper does not explicitly discuss this issue. 
It is however clear that because of the decomposition = Cii^'^" + C'^'*^°) the loop corrections 
would have the form 

with IR cutoff go ' 

(45) 

and thus the total correlator is: 

A'in, r?)tree +A'{k., 7?)loop, with IR cutoff go • (46) 

It is straightforward to see then, that, although the "IR-fields" are resummed in k, equation 

Actually, |20l 121] discuss the spectrum, which is the two-point function rescaled by fc^ , but this is not relevant 
here. 
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()35p is indeed just given by a coordinate transformation and the final outcome amounts in 
effect just to the use of an infrared cutoff qq in the loops. That is, there is in effect no difference 
with the usual "box cutoff" approach, with all the issues that we have discussed before. 

It follows that there are big differences with our current proposal: we do not introduce any 
infrared cutoff whatsoever, not even an "observational" one like qo, nor we decompose the fields 
as C = Cir'^'^" + C''^'^"- What we do is really to define a new perturbation theory on the basis of 
the quantum expectation value h = {e^^e"') of the full quantum operator to define scales, instead 
than the standard perturbation theory around the background metric. The new perturbative 
correlators are given by the series Eq. ()4ip . and the higher order corrections to correlators 
are non-trivial and do not reduce to the use of a cutoff in the old perturbation theory. It is 
straightforward to generalize to all iV-point functions for both Q and 7. 



4.4 Example at one-loop 

We now give a detailed example at one-loop that our procedure do cancel all IR divergences (Log- 
cutoff enhancements) while redefining the perturbative series, in the theory with the Lagrangian 
in ([22]) . To keep the example simple, let us consider a pure de Sitter case, so that C = 0,(7 = 
a{x,t). In this case, the only physical metric perturbations are gravitons. We wish to compute 
the IR-safe two-point correlator {a{k,ri)a{k ,rf)), at one loop, following our recipe culminating 
in equations ([^T]). 

We start by computing the one-loop correlator using the standard in-in perturbation theory. 
We recall here the known results from Ref . [llj . The one-loop in-in loop corrections are given by 
the diagrams in Figs. [T]and[5J The relevant interaction Lagrangians are obtained by expanding 
the full Lagrangian in perturbation up to second order, and they read 

= ^Y^d,adja , (47) 

C^ = -^Y'j^'diadja, (48) 

where a is scale factor of the de Sitter background. The computation is performed by using 
Eq. (|12p . and the details are quite complicated and of little interest here. We will give some of 
the intermediate steps, details can be found in Refs. [9l [11]. From the diagram in Fig. [U we 
obtain two contributions 

d^q p dm dm .-^,J,J,fc,JT 



A,ir^) = -8Re / ^ J ^JW^ j ^J^2^^,Mmk' k'k'w,{m.m)Wt-,{m,m) 

X VFfc (7?, r?i)VFfc (77,772) (49) 

k {Hriif 

(50) 

where Wa{rj,r]') is given by Eq. From Fig. [2l we obtain 



5fc(r?) = 4Re j ^ ]^ ,;(g-;)fc^A;^fc'=fc'M^,( 



(51) 



Ck{il) = Re 

The definition of the sum over graviton polarizations uJij^ki is 

^ij,kii?) = ^4j{^^ki{(l) = ^ik^ji + ^ii^jk - ^ijhi (52) 

s 

+ Sijqkqi + himj - SikqjQi - ^iiQjqk - Sjkmi - ^jimk + miWi 



^In fact, equation H44|) is similar to (|43p . that is difFerent from our (|4T 
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where q is the unit vector in the direction q, and Cij is the graviton polarization tensor. After 
evaluating the integrals and regularizing it with IR and ultraviolet cutoffs, we get 

{a {k,r])a{k', 7])) one-loop = Ak{r]) + Bk{r]) + Cfc(r?) 

2 . . . .H'^ 1 



+ {finite terms and UV corrections} + O ( j (53) 



Therefore, adding the three-level result 



{a{k,v)a{k',7l))tree = {27T^)6{k + k') ^^(1 + k\'') , (54) 



the total result for the two-point correlator is 



16. 



A{k,vi) = {a{k,rj)a{k, r]))trcc + {cr{k,r]) a {k,rj)) one-loop 

{finite terms and UV corrections}] -|- O (^'j^^ 



(55) 



We have reinstated the Planck mass to show explicitly the suppression factors H'^/Mp. Note 
that |7|2 ~ HyMl 

We now wish to apply our procedure to obtain the IR-safe correlator up to one-loop order, 
starting from this result. We follow the steps in section 14.31 The first one is to compute the 
terms of the series in Eq. ()4ip up to the order of interest. Since we are working up to one-loop, 
this means up to the order ~ which is quadratic in 7. Looking at Eq. (f55|) . we see that we 

need to compute up to the second order derivative of {a{k,r])a{k,7]))tj.(,e, and no derivatives of 
{a{k,r])a{k,r])) 

one-loop • The second order expansion in 7 gives 
^{l'^l^^)d^rjd^u{(T{k,vMk',v))trec\^2=-^2 = ^ ^{l%^%^jdk2{a{k,ri)a{k',r{))treo\^2=-^2 

id 

+ \ Y.^l"^^^^)h^j^khdl2{a{k,ri)a{k',r]^ 

(56) 

After some manipulation, we obtain the contributions 

:jy'(7''7z^>^ifcj42(0-(^,'?)f^(fc',??)>tree|;.2=fc2 = - o^Og{Al r)—^ -=3^ r]^ + 

+ {disconnected finite terms and UV corrections} (57) 



lY,{Y^^''^)kikjkkkidl2{a{k,r])aik',r]))tree\f^2^-^2 = ^log(Ajfl)^ \_s k\'^+ 

+ {other disconnected finite terms and UV corrections} . (58) 
'^^In [11] . only the result for kr] ^ were presented. 
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We now apply Eqs. (fH]) . subtracting Eqs. ([571 158|) form Eq. ([55]) . to obtain 



^ _ H 1 —2 

Am.sMk,k' ,ri) = {2-jr^)6{k + k') — j — + ??^)+ 




) 



(59) 



As we can see, the dependence of log(A/jj) has disappeared and there is no trace of any other IR 
regulator. The new finite terms and ultraviolet corrections are given by the difference between 
those in equation ([55]) and the disconnected contributions in (fST]) . ([58]) : they are non-zero as 
is not given by disconnected diagrams. 

4.5 A^-point functions and non-Gaussianities 

What we have said about the two-point function can be extended to the most general case 
of A^-point functions. It is interesting then to discuss the loop correction to the three-point 
function and the bispectrum, which is the first indication of non-Gaussianities in the primordial 
perturbations. 

The latter ones are usually quantified with a set of parameters. The one which is related to 
the three-point function is called Jnl, having defined 



The parameter fi\fL will receive corrections at one- loop. In particular, in the standard approach, 
it suffers from IR Logarithmic-enhancements (although mitigated by some dependence on the 
spectral index 1 — re) [11]. Such corrections can make non-Gaussianities much more pronounced 
theoretically than what is actually measured, see 0. 

Instead, if we use our IR-safe proposal to calculate the three-point function corrected up to 
one loop, one does not find such Logarithmic-enhancements, but only the ultraviolet corrections 
(which are renormalized). We will not compute them here as it goes beyond the scope of this 
paper. 

5 Conclusion 

In this work we have proposed a method to solve the IR issues of cosmological correlators. The 
solution we propose here is based on the physical principle - every observable (correlator) should 
be defined in terms of a local measurement. Of course, when one in practice approximates the 
ensemble averaging of the correlators with the spatial averaging, one introduces by definition a 
certain degree of non-locality. We have discussed the connection between the IR issue and this 
approximation at the end of section [3l 

Our proposal consists in a modification of the formalism enabling us to ask physical questions. 
In particular, we emphasize that quantities must be defined locally by local observers and we 
implement this principle by defining a new perturbation theory in the in-in formalism. 



We do not see any breaking of perturbation theory due to IR corrections [32] , as our definition 

of the perturbation theory is IR-safe. Furthermore, there is no ambiguity related to the choice 
of regularization of the IR divergent integrals, as the final result does not depend on it. 



(60) 



and written the bispectrum as 




(61) 
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